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Abstract 

Based on the analogy with the quantum mechanics of a particle propagating 
in a complex potential, we develop a field-theoretical description of the sta- 
tistical properties of a self-avoiding polymer chain in a random environment. 
We show that the account of the non-Hermiticity of the quantum Hamiltonian 
results in a qualitatively different structure of the effective action, compared 
to previous studies. Applying the renormalisation group analysis, we find a 
transition between the weak-disorder regime, where the quenched random- 
ness is irrelevant, and the strong-disorder regime, where the polymer chain 
collapses. However, the fact that the renormalised interaction constants and 
the chiral symmetry breaking regularisation parameter flow towards strong 
coupling raises questions about the applicability of the perturbative analysis. 
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I. INTRODUCTION 



The problem of a polymer in a random environment is among the most interesting in 
statistical physics. It has been known for a long time that the mean square end-to-end 
distance of a pure self-avoiding walk (SAW) of length L obeys the scaling law (r 2 ) ~ L 2u , 
where v m 0.59 in three dimensions (for a Gaussian random walk, one has the classical 
exact result v = 0.5 in all dimensions). The question of how this scaling behaviour is 
affected by external impurities has attracted considerable research effort for more than 
a decade H5H17H- In his pioneering work, Harris || argued that, treated perturbatively, 
quenched disorder is irrelevant, and, therefore, no modification of the critical exponent v 
should be expected (see also Ref. f7]). This conclusion found support in the Monte-Carlo 
simulations on weakly diluted lattices ||. The opposite case of strong disorder has also 
been studied, both numerically || and analytically flJT0||. Edwards and Muthukumar || 
and Cates and Ball |TTJ predicted that a Gaussian chain placed in the field of impurities 
would collapse to a localised state, in which (r 2 ) — > const at L —>■ oo. The localisation 
breaks down if one introduces a weak repulsive interaction between the monomers in the 
chain. In this case, it was found ||TT| that the polymer behaves as a free random walk with 
v = 0.5. Later, a crossover between the regimes of weak (y ph 0.59) and strong [y = 0.5) 
disorder was predicted to occur at some critical concentration of impurities [O. However, 



these results seem to contradict the conclusions of a number of other authors [|14yi5|l , who 
argued that, at any concentration of impurities, the scaling of a long chain with excluded 
volume interactions is controlled by the critical exponent of a directed random walk, v = 2/3. 
Overall, a comprehensive theory that would describe the effects of disorder on self-avoiding 
polymers is still missing. 

Previous studies of the interplay of disorder and excluded volume interactions have been 
largely based on qualitative arguments, although a number of mathematical techniques 
(including the variational methods [ 11], 12]] , replica field theories 0,|13|,|14|] and the real-space 
renormalisation group [|16j) have also been used. A more general framework would therefore 
be beneficial to a better understanding of these effects. We believe that such a framework 
could be designed by employing a field-theoretical approach, which draws on the connection 
between the statistical mechanics of a polymer and the quantum mechanics of a particle 
in a complex random potential. Despite its apparent simplicity, the use of this method has 
been hampered by the fact that, as the Hamiltonian of the particle in a complex potential is 
non- Hermit ian, it is impossible to represent the Green function in the form of a convergent 
functional integral (see below). A similar problem has been encountered, and successfully 
resolved, in the spectral theory of non-Hermitian operators. The latter has attracted great 
interest in recent years. A variety of applications have been identified including the study 
of anomalous diffusion in random media |18|, scattering in open quantum systems ||19|| , 



neural networks [2~0|, and the statistical mechanics of flux lines in superconductors |2T]. The 



problem of a self-avoiding walk without impurities, which can be mapped onto the quantum 
mechanics of a particle propagating in a random imaginary potential, has also been analysed 
in this context |22| . These studies have led to the development of a new technique, based on 
the representation of the spectral properties of non-Hermitian operators through an auxiliary 
Hermitian operator of twice the dimension 2H|| , which serves as a good starting point for a 
field-theoretical approach. 
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The main purpose of the present paper is to derive a consistent field-theoretical formu- 
lation of the problem of a self-avoiding polymer chain in a random white-noise potential, 
using the methods of non-Hermitian quantum mechanics, which is done in Section PJ. The 
large-scale behaviour of this model is studied perturbatively in Section OJ by means of the 
momentum-shell renormalisation group (RG) in D = 4 — e dimensions. Section concludes 
with a discussion of the results obtained and of possible limitations on the validity of the 
perturbative approach. 



II. DERIVATION OF THE FIELD THEORY 

Let us consider a continuum self-avoiding chain of length L, with one end fixed at the 
origin. Then the probability to find the other end at a point r can be expressed as a path 
integral (Edwards model p4] , p5| ): 

<-*M=* _ r 1 t L (d*(s)\ 2 r L 



P(r,L)= f r Dx( S ) expj-^ f ds ^ - / ds V^s)) 
A(o)=o I 2a Jo \ ds J Jo 

-y [ d Sl | o L ^ 2 5(x( Sl )-x(s 2 ))}. (1) 



The first term in the exponent corresponds to the entropic contribution (the length a, called 
the Kuhn length, is a microscopic parameter with the physical meaning of the monomer size, 
which provides a natural ultraviolet cutoff scale). The second term is the potential energy 
of the chain in an external potential, which is assumed to be a Gaussian distributed random 
function with the correlator Vi(r)l / L (r / ) = 7j5(r — r'). The last term takes into account 
the excluded volume effects (72 > 0), the limit -y 2 — ► 00 describing the situation where the 
intersections of different fragments of the chain are penalized by an infinite energy barrier 
and thus completely forbidden. This interaction term can be decoupled by introducing an 
auxiliary white-noise potential. Using the identity 



eX F>{~-^/ dsi J ds 2 5(x(si) - x(s 2 ))| = ^exp | — dsV 2 



:*oo) 



where the angular brackets denote averaging over a Gaussian distributed random field V^r) 
with correlator (V^r)^!"')) = 7 2 5(r — r'), we can represent Eq. (JJ) as an averaged Feynman 
propagator of a fictitious quantum particle moving in a complex random potential: P(r, L) = 
(U(r,L)), where U(L) = exp(-LH) with the Hamiltonian 

H = -V 2 + Vi(r) +iV 2 {r) (2) 

(Vi s 2 (r) are independent Gaussian random fields). Note that H is non-Hermitian and may 
therefore have complex eigenvalues. 

It is convenient to use the "energy representation" , which is achieved by introducing the 
Green operator ClS db function of z = x + iy, 

= = £!**>— (3) 

z — ri u z — z k 
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where \Rk) and (L&| are the right and left eigenfunctions of H, and z^ denote the complex 
eigenvalues. Using the identity dz~ x jdz* = n5 2 (z) = ir8(x)5(y), one can relate U(L) to 
g(z), and express the end-to-end probability distribution of a self-avoiding chain in a given 
distribution of impurities in the form 

P{r,L\V x ) = lJd 2 z exp(-zL)£-(g{r,z)), (4) 

where the integration runs over the entire complex plane. It should be noted that the 
standard field-theoretical methods, based on the representation of an inverse matrix in the 
form of a Gaussian functional integral, cannot be directly used for the calculation of g(z), 
the reason being that, as H may have eigenvalues anywhere in the complex plane, it is 
impossible to guarantee the convergence of the functional integral. As pointed out in Ref. 
|26fl , a naive attempt to calculate the functional integral by the analytical continuation from 



the region where it converges to the whole complex plane fails. The problems are revealed 
by representing the density of complex eigenvalues through the identity 

1 d 

p(z)=Y, 52 ( z - z k) = --^^V9(z), (5) 

k 71 oz 

wherein the Green function is shown to be non-analytic everywhere in which the density of 
states is non-vanishing. 

To circumvent these difficulties, a representation has been introduced , in which 

the complex Green function g(z) is expressed through an auxiliary Hermitian operator, which 
in our case has the form 

G-\z) = ( ^ _° £ f z ~ H ) = (x + V 2 - \4(r)) a.-iy- V 2 (r))a 2 , (6) 

where cij are the Pauli matrices. A relationship between g and G is straightforward: 

g(z) = G 21 (z). (7) 

Using the replica trick, the matrix Green function G(z) can be written as a functional 
integral over 2n-component complex Bose fields (in the limit n — > 0): 

G*(r, z) = -^limj V 2 cp a exp {z J d D r ^ a {G~\z) + ^o)v a } ^(r)^f(O) (8) 

(o"o is the unit matrix, summation over repeated replica indices is assumed), where 



9/ \ ' ~a II 11 



Due to the Hermiticity of and the presence of the term with 77 ("regulator"), the 

functional integral is well defined and convergent. Note that, although the matrix Green 
function @ possesses a chiral symmetry, a^G~ Y {z)a^ = —G^ 1 (z), this symmetry is broken 
in Eq. (g) if 77 ^ 0. 



4 



Eqs. (0) and (§) allow one to average g(z) over V2 and calculate the polymer partition 
function in a given configuration of the external disorder. In order to obtain physically 
observable quantities, one has to average the free energy F(r,L) = lnP(r, L) over the 
quenched random field V\, which can be done using the replica trick once again: 



F(r, L) = lim 



P m (r,L) - 1 



m-*0 m 



1 f r m d 2 z ( rn \ m d " l I 

iim n - / n — ex p - L ^>« ii^n^a-i • (9) 



It is thus necessary to introduce a second set of replica indices and to integrate over 2nm- 
component fields ^ l aa ij) (i = 1,2; a — 1, ...,n; a — 1, ...,m). Averaging over Vi, we obtain 
from Eqs. (0) and (g): 

n^Xj) = / ^V a>a e^UV-J n^LW^ ; ;(0) (n,m - 0), (10) 
where the effective action has the following form: 

iS[(f] = J d D r {i{v\ a VlV 2 Va,a) + ix a(<Pl, a VlV a , a ) ~ Waivla^V^a) - V^cxW a,a) 
Tl i t T2 I t 



y(<a^l'^a,a)(<a 0r l^,a)}- (H) 

This action is different from the replica n-vector model, used previously in the RG analysis of 
the problem of a SAW in a random medium ||7|. [T3| . [l4| , because of the double dimensionality 
of the fields involved and a larger number of the coupling constants (three instead of two). 
The additional term proportional to 7^ cannot be obtained in a formal derivation of Eq. 
(pT|), but should be added to the long- wavelength effective action for consistency, as we 
shall see below. Note also that we have kept the term with 77 in Eq. ( |TTD (the importance 
of this term will become clear shortly). The asymmetry between the ways the Latin and 
Greek replica indices appear in Eq. (|TT| ) reflects the differences in the nature of the random 
potentials : V\ describes the external quenched disorder, while V2 is a fictitious annealed 
random field. 



III. RENORMALISATION GROUP ANALYSIS 

The long-distance properties of the field theory ( [TT]) can be investigated using the mo- 
mentum shell RG approach and the e-expansion near the upper critical dimension D c = 4. 
Using the standard procedure, consisting of the separation of "slow" and "fast" degrees of 
freedom followed by a rescaling of lengths and fields |f29fl , we obtain the flow equations in 
one-loop order: 
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dlnxn 1 , , s 

" 2 + —( 7l -7 2 + 70, 



d£ 8vr 

d^Va 1 , ,* 

2 + ^wi-72 + 7i), 



d£ 8tt 2 
cilnr? 1 . 

V = 2 + ^ (7l + 72 + 7l) ' 

^ = 67i + ^( 2 7? - 7i72 + 7i7i), (12) 

^72 1 / 2 /\ 

^- = ^72 + ^(3 7 i 7 2 - 7 2 + 3 7 27i), 

= e7 i + 4^ 3 ^ + ^ - 7271 + 2 7 f), 

where £ = ln(L/a) is the logarithmic RG parameter. The real and imaginary components of 
the complex "energy" z = x + iy are renormalised in the same way. It also follows from Eqs. 
(|12|) that, even though initially 7^ = 0, it acquires a non-zero value under renormalisation. 
Introducing a different set of independent variables, 

u = 71 -72 + 7i ; ^ = 72, w = i Xl 

we rewrite Eqs. (|l|) as 



dhix 
d^ 


1 d In y 1 

= 2+ 8^' di =2+ 8n^ 


(13) 


din T] 
d£ 




(14) 


du 
*T 


1 2 

eM + 0^2 M ' 


(15) 


dv 

d£ = 


ev + ^(3uv + 2v 2 ), 


(16) 


dw 
~dj ~ 


ew H ^(3uu; + u 2 + 2t>io — w; 2 ). 

47T 2 


(17) 



Note that the flow equations for x, y and u are decoupled from those for v and w. 

Let us neglect the renormalisation of rj, i.e. put 770 = (the consequences of allowing for 
r)o 7^ will be discussed below). As seen from Eq. (H), the main contribution to the average 
free energy at L — > 00 comes from x,y ~ 1/L, which provide the initial conditions for the 
equations (0). The scale R c , at which the renormalised values of x and y become of the 
order of unity, represents the correlation length of our field theory and should be identified 
with the average size of a polymer of length L. Integrating Eqs. (|i"3|) with respect to £ from 
£ = to £ c = \n(R c /a), we obtain the equation which relates the polymer size to L: 

>■+=/*> ( 2 + |f). (18) 
a Jo \ 8ir z J 

The behaviour of essentially depends on the bare values of the coupling constants 7^2 
(the bare value of j[ is zero). One should distinguish between the two possibilities: In 
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the weak disorder regime (71 < 72, i.e. uo < 0), u(£) flows towards the stable fixed point 
u* = — 2ir 2 e. After substitution in Eq. fll8|) , we obtain R C (L) ~ L v with z/ « 1/2 + e/16. 



This exponent coincides with the result for a pure SAW |lj which means that weak disorder 
is irrelevant for the large-scale properties of polymers. This result is in agreement with Refs. 

BHD- 

In the strong disorder regime (71 > 72, i.e. uo > 0), the solution of Eq. (|15|) 
has a pole at £ = £*, where 

r = -i-ri+— V (20) 

Substituting tt(£) ~ (£* — £) 1 i n Eq. Ql5[), we find that, at L — > 00, R C (L) — >■ ae^*. Thus, 
the polymer size tends to some constant independent of its length. Although this conclusion 
coincides with the results of Refs. P JT^JT7[ ], the applicability of our one- loop RG equations 
in this strong-coupling regime is limited. 



IV. CONCLUSIONS AND DISCUSSION 

Our results show that a proper account of the non-Hermiticity leads to significant changes 
in the structure of the effective field theory for self-avoiding walks in random media. The 
application of the one-loop renormalisation group analysis to the action with the doubled 
number of the degrees of freedom allowed us to predict the existence of a phase transition 
between the weak disorder regime (in which the scaling behaviour of the polymer is not 
affected by the external randomness), and the strong disorder regime (in which the polymer 
collapses into a localised state). This conclusion seems to be in agreement with one of the 
scenarios considered previously, notably in Refs. [f7|,|16 . 



However, our RG calculations also indicate that there might be some hidden flaws within 
the perturbative approach, one of which is revealed by taking into consideration the infinites- 
imal regulator 7/, which explicitly breaks the chiral symmetry of Eq. (||). It follows from Eq. 
fllED that, in the weak disorder case, i>(£) has a pole at £ = £*, where 



C = - In 1 + 2vr Vr° " ' . (21) 




Av - 


\u \ 


(2^o- 


kl) 2 



After substitution in Eq. fll4|) and integration over £, we find that r/(£) grows under renor- 
malisation, 7/(£) ~ 7/0(1 — £/C) -< ^ 2 5 an d diverges at some scale R C)V ~ ae&. Although 7/ 
does not appear explicitly in the flow equations for x, y and u, its singular RG behaviour 
may be an evidence of the spontaneous breaking of the chiral symmetry |I0|. In the field- 



theoretical language, it would manifest itself as the appearance of the anomalous pairings 
such as (ip l '*Lp l ) and (ip 2 '*ip 2 ), which could change the long-distance behaviour of our sys- 
tem at scales larger than R CjV . It might be the case that a seemingly trivial technical trick 
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of doubling the degrees of freedom actually reveals deep physics related to the presence 
of the additional continuous symmetries in the system, whose effects would otherwise be 
completely lost. 

Furthermore, despite the presence of a stable fixed point in the RG equation fll5|) for u 
in the weak disorder regime, the coupling constants 71, 72 and 7^ are all renormalised to 
infinity at large scales (see Eqs. ( |16D and (|17D ). The existence of the underlying RG flows 
towards a strong coupling regime can cast doubts in the validity of the scaling description 
itself (the authors of Refs. [|l4] , |l"5f suggested that the behaviour of the system could in- 
stead be determined by a strong-coupling fixed point, which is inaccessible by means of the 
perturbation theory). It is interesting that an analogous situation takes place even in the 
absence of an external disorder (i.e. at 71 = 0), in which case 72 and j[ both go to infinity, 
but 72 — 7i — > 27r 2 e. One of the possible explanations is that an important role could be 
played by the so-called "tail states" in the spectrum of the Hamiltonian (|j), which would 
significantly change the observable large-scale behaviour of the system. The properties and 
physical manifestations of the "tail states" of Hermitian operators have been extensively 



studied since the pioneering works of I. M. Lifshitz pi| . A generalisation of this theory to 
the case of non-Hermitian random operators presents an interesting theoretical challenge 
with numerous potential applications. 

Finally, we would like to emphasize that, since at every step of the calculations we deal 
with convergent functional integrals, our field-theoretical model of SAWs in random media 
is not only capable of reproducing some of the previous results on the perturbative level, but 
also seems to suit well for the study of non-perturbative effects, such as the chiral symmetry 
breaking and/or the contribution of the spectral "tails". We leave the investigation of such 
effects for future work. 
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